
caixeta <- read.csv("caixeta.csv", as.is=T)
names(caixeta)
table(caixeta$especie)



hist(caixeta$cap/pi) 


agr=aggregate(caixeta$h, by=list(local=caixeta$local,arvores=caixeta$arvore), FUN=mean) ##rvores com mais de um fuste so agregadas pela altura mdia, mantendo-se a informao do local.

chauas=agr[grep("chauas", agr$local),]

jureia=agr[grep("jureia", agr$local),] 

retiro=agr[grep("retiro", agr$local),] 

par(mfrow=c(3,1)) #Visualizao dos trs histogramas na mesma pgina

hist(chauas$x) #Frequncia de rvores por altura em chauas

hist(jureia$x) #Frequncia de rvores por altura em jureia

hist(retiro$x) #Frequncia de rvores por altura em retiro

par(mfrow=c(1,1)) #Visualizao de grficos para um por pgina.


#H diferenas entre as estruturas (distribuio de tamanhos) dos caixetais?

summary(chauas)
summary(jureia)
summary(retiro)


par(mfrow=c(3,1)) ##Ajusta a visualizao dos trs histogramas na mesma pgina.

hist(chauas$x, ylim=c(0,0.015), xlim=c(0,500), prob=T, breaks=21) ##breaks  necessrio para igualar a diviso de barras dos outros grficos.

lines(density(chauas$x), col="blue") ##Ajuste de densidade probabilstica.

curve(expr = dnorm(x,mean=mean(chauas$x), sd=sd(chauas$x)), add=T, col="red") ##Curva normal.

hist(jureia$x, ylim=c(0,0.015), xlim=c(0,500), prob=T) #Ajuste o eixo x com o eixo do grfico anterior.

lines(density(jureia$x), col="blue") ##Ajuste de densidade probabilstica.

curve(expr = dnorm(x,mean=mean(jureia$x), sd=sd(jureia$x)), add=T, col="red") ##Curva normal.

hist(retiro$x, ylim=c(0,0.015), xlim=c(0,500), prob=T) #Ajuste o eixo x com o eixo do grfico de chauas.

lines(density(retiro$x), col="blue") ##Ajuste de densidade probabilstica.

curve(expr = dnorm(x,mean=mean(retiro$x), sd=sd(retiro$x)), add=T, col="red") ##Curva normal.

par(mfrow=c(1,1)) ##Ajusta a visualizao dos trs histogramas na mesma pgina.
##Apesar de apresentarem mdias diferentes e medianas no discrepantes em relao s primeiras, os locais avaliados no apresentaram grandes alteraes na estrutura.
